Positioning systems in Minkowski space-time: Bifurcation problem 

and observational data 



Bartolome Collj^j Joan Josep Ferrandoj^j and Juan Antonio Morales-LladosE^] 

Departament d'Astronomia i Astrofisica, Universitat de Valencia, 46100 Burjassot, Valencia, Spain. 

In the framework of relativistic positioning systems in Minkowski space-time, the determination of 
the inertial coordinates of a user involves the bifurcation problem (which is the indeterminate location 
of a pair of different events receiving the same emission coordinates). To solve it, in addition to the 
user emission coordinates and the emitter positions in inertial coordinates, it may happen that the 
user needs to know independently the orientation of its emission coordinates. Assuming that the 
user may observe the relative positions of the four emitters on its celestial sphere, an observational 
rule to determine this orientation is presented. The bifurcation problem is thus solved by applying 
this observational rule, and consequently, all of the parameters in the general expression of the 
coordinate transformation from emission coordinates to inertial ones may be computed from the 
data received by the user of the relativistic positioning system. 



I. INTRODUCTION 

To locate the usertQ of a Global Navigation Satellite 
System (GNSS), several geometric methods and algebraic 
algorithms have been developed in the past [IH5] that are 
still in use [51 [7]. Basically, the algebraic statement of 
the location problem is rather simple: to find the events 
where the emission light cones of four broadcast signals 
intersect. Of course, this idea is implicit in the Bancroft 
algorithm [5] and other similar ones [3]. In fact, Abel 
and Chaffee |H [5] used Minkowskian algebra to state 
the problem properly, making apparent that the more 
Lorentzian a description is, the more clear algorithm is 
performed. 

However, in a full relativistic framework (cf. [81414] ) . 
and even in the case of the fiat space-time, an explicit 
form of the solution of the location problem for arbitrary 
emitters has not been obtained until recently [15] . 

In |15j , an exact relativistic formula giving the inertial 
coordinates of an event in terms of the received emission 
coordinates is obtained. This formula applies in all the 
emission coordinate region and involves the orientation of 
the emission coordinates of the user. Nevertheless, there 
exists an inherent limitation on the applicability of this 
formula: only the users in a certain region (named the 
central region, see Sec. IV B I of a positioning system can 



obtain the orientation from the sole standard emission 
data, that is to say, from the sole set of the positions of 
the four emitters in inertial coordinates and of the emis- 
sion coordinates of the user. Consequently, only these 
restricted users are able to locate themselves in inertial 
coordinates. 



Here, assuming that the users out of the central region 
may observe the relative positions of the four emitters 
on their celestial sphere, we will give a simple rule allow- 
ing any user of the positioning system to locate itself in 
inertial coordinates. To show that, we will see that the 
orientation of the emission coordinates of a user is related 
to the relative positions of the emitters of the positioning 
system on the celestial sphere of the user. 

In building current GNSS models, the usual assump- 
tion consists in picking out an approximate numerical 
solution. But, because gravitational effects are not taken 
into consideration at the considered leading order, one 
should start from the best accurate solution that nowa- 
days we know. Such a solution is precisely the simple, 
exact, and covariant formula, found in Ref. |15j and im- 
proved here, giving the location of a user of a relativistic 
positioning system in Minkowski space-time. 

Let us remark that our result not only concerns GNSS 
around the Earth, but also general (relativistic) position- 
ing systems anywhere in the Solar System or elsewhere. 
It is true that for most, but not all, of the present ap- 
plications of the GNSS the users are near the Earth's 
surface. Therefore, they are usually in the central region 
of the satellites they detect, so that additional data (and 
in particular our observational rule) are not necessary. 
But for other applications of the GNSS as well as for 
general positioning systems, our observational rule may 
be a simple way to solve the bifurcation problem and 
hence, the location one. 



A. Outline of the paper 
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The word "user" here denotes any person or device able to re- 
ceive the pertinent emitted data from the relativistic positioning 
system and to extract from it the corresponding information. For 
short, we shall refer to the user as "it". 



The paper is organized as follows. In Sec. [ITJ the in- 
ertial coordinates of a user are expressed in terms of the 
emitter configuration and the orientation of the position- 
ing system. This provides a covariant formula for the 
transformation from emission to inertial coordinates. An 
analysis of the solution in terms of the configuration of 
the emitters is also presented. In Sec. |III[ some proper- 
ties of the border between the two emission coordinate 
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domains are obtained, and an observational rule to de- 
tect it is remembered. Section Hvl is devoted to define the 
genuine regions and coordinate domains involved in the 
problem. We stress the geometrical meaning of the coor- 
dinate transformation formula in connection with these 
regions. In particular, we show that in the central region 
of the positioning system the orientation is computable 
from the sole standard emission data. In Sec. |V]we dis- 
cuss the bifurcation problem (nonuniqueness of solutions 
in the determination of the location) which is related to 
the existence of regions whose events can not be located 
from the sole standard emission data. We give an obser- 
vational rule to solve the above indetermination problem. 
This rule allows us to determine, at any event in the emis- 
sion region, the orientation of the emission coordinates 
of the user from the observational data of the relative 
positions of the emitters on the celestial sphere of any 
user at this event. The concluding Sec. |VI| is devoted to 
summarize and discuss the results. The used notation is 
explained in an appendix. 

Some preliminary results of this work were presented 
at the Spanish Relativity meeting ERE-2010 [T6] . 

B. Relativistic positioning terminology: Brief 
compendium 

As pointed out in Ref. [5], relativistic positioning sys- 
tems [8 13J and the emission coordinates |14[ I17j they 
realize are essential elements to develop the relativistic 
theory of the GNSS. Starting from scratch, we present 
here a compendium of basic definitions about this specific 
subject. Anyway, we consider these definitions necessary 
not only to make this paper self contained, but also as 
an incipient piece of concepts to deal with GNSS in a full 
relativistic perspective. 

Relativistic positioning system: set of four emitters A 
(A = 1,2,3,4), of worldlines 7_a(t" 4 ), broadcasting their 
respective proper times t a by means of electromagnetic 
signals 

Emission coordinates of an event: the four times {t a } 
which are received at each event reached by the emitted 
signals]^] 

Configuration of the emitters for an event x: set of four 
events {7^ {t a )} of the emitters at the emission times 
{t a } received at x. 

Emission region: set TZ of events reached by the four 
signals broadcast by the positioning system. Every x G TZ 
is labelled with the corresponding emission coordinates 
{r A }. 



2 For simplicity the proper time is taken here, but any other time is 
valid. For example, the Global Positioning System (GPS) broad- 
casts the GPS time, a time which, roughly speaking, coincides 
up to a fixed shift with the International Atomic Time (TAI), a 
sort of mean proper time on the Earth surface. 

3 Emission coordinates have received different appellations in the 
past (see |18| for a brief and critical account). 



Characteristic emission function: map that to every 
x € TZ associates its emission coordinates, that is Q(x) = 
(r A ). 

The characteristic emission function describes the ac- 
tion of a positioning system and, hence, represents it. 

Emission coordinate region: subset C of the emission 
region TZ where the gradients dr A are well defined and 
linearly independent. 

The emitter worldlines arc excluded from C because 
every dr A is not defined at the emission event ja {t A ) 
(this event being the vertex of the emission light cone 
t a = Constant). 

Orientation of a relativistic positioning system at the 
event x: orientation of its emission coordinates at x. It 
is given by the sign e of the Jacobian determinant jq (%) 
of at x, e = sgnjQ^x). 

In terms of the gradients of the emission coordinates, 
one has 

e= sgn^idr 1 A dr 2 A dr 3 A dr 4 )} (1) 

where * stands for the Hodge dual operator, and A is the 
exterior product (see Appendix [A] for transcription into 
index notation). 



II. THE LOCATION PROBLEM IN 
MINKOWSKI SPACE-TIME 

Suppose a given specific coordinate system {x a } cover- 
ing the emission region TZ, let 7^ {t a ) be the worldlines of 
the emitters referred to this particular coordinate system, 
and let {t a } be the values of the emission coordinates re- 
ceived by a user. The data set E = {ja{t a ),{t a }} is 
called the standard emission data set. 

The location problem with respect to E, also called the 
standard location problem for short, is the problem of 
finding the coordinates {x a } of the user from the sole 
data E. 

In Ref. [T5] , the above standard location problem was 
analyzed for arbitrary relativistic positioning systems in 
Minkowski space-time, assuming that the specific coordi- 
nate system {x a } is an inertial one. There, the explicit 
expression x a — K a (r A ) was found, giving the coordinate 
transformation from emission coordinates to inertial ones 
(Eq. {3} below). 

Particular simple cases have already been studied: con- 
sidering a 2-dimensional [TT?rl2"T] or a 3-dimensional [55] 
space-time, or for special motions of the emitters in the 
Schwarzschild geometry from analytical |23j and numer- 
ical |24j approaches. For a recent approach to emission 
coordinates using the integration of the eikonal equation, 
and some numerical simulations, see also Ref. |25j . 

In this and the following sections, we are mainly deal- 
ing with relativistic positioning systems in Minkowski 
space-time. 



A. Covariant expression of the solution 

From now on, we shall suppose that any user in the 
emission coordinate region C receives the standard emis- 
sion data set E. Let us denote by x the position vector 
(with respect the origin O of this inertial system) of an 
event P in the emission region 1Z, x = OP. If a user 
at P receives the broadcast times {t a }, ja denote the 
position vectors of the emitters at the emission times, 
1A ee 1a {t a ). Then 



m A = x - ja, 



(A=l,...,4), 



(2) 



are future-oriented light-like vectors that represent the 
trajectories followed by the electromagnetic signals from 
the emitters ja {t A ) to the reception event x G 1Z (see 
Fig.0). 

In the standard emission data set E, the emission data 
{t a } received at x are the emission coordinates of the 
event ieK and were broadcast when the emitters were 
at the events {"/a(t a )}, the configuration of the emitters 
for the event x. Generically, these four events determine 
the configuration hyperplane for x^ 

For the events x in the emission coordinate region C, 
the transformation x = k(t a ) from emission to iner- 
tial coordinates is locally well defined. In [15], we have 
obtained a covariant expression of this transformation, 
given by the following formula: 



74 + y* - 



vlx 



(y* ■ x) + £\/(y* ■ x) 2 



ylx 2 



(3) 




FIG. 1. The emission of an electromagnetic signal from a 
satellite 7a (t A ) at proper time t a , and its reception by a 
user at x = OP. These events define the future pointing null 
vector m A - 



4 Here, we always consider that the emitter configuration is regu- 
lar, i. e., the four events {7a ( rA )} are noncoplanar. Nonregular 
or degenerate configurations are considered elsewhere (for some 
remarks concerning these situations, see Refs. [U0). 
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FIG. 2. If we choose the emitter four as origin (reference 
emitter), the relative positions of the others (referred emit- 
ters) are e a =7a — 74, a — 1, 2, 3, and the position vector of 
the event P is 7714. 



where 74 (t 4 ) has been chosen as the reference emitter^] 
Quantity is given by 



1 



y* 



(4) 



where \ is the configuration vector 

X = *(e x A e 2 A e 3 ) (5) 
and H is the configuration bivector 

H = e 2 A e 3 + Q 2 e 3 A e\ + Q 3 ei A e 2 ) (6) 
with (see Fig. ^ 

^a = ^(e a ) 2 , e a =7 a -7 4 , (a =1,2, 3), (7) 

and where £ is any vector transversal to the configuration, 
£ • X 7^ 0, and i{£,)H stands for the tensor contraction of 
£ and the first slot of H (see Appendix [A} . 

Quantity e is the orientation of the positioning system 
at x, that is now equivalently expressed as 



sgn[*(mi A m 2 A m 3 A 7714)]. 



(8) 



It is worthy to remark that x an d H are determined by 
the relative positions e a = 7 a — 74 associated with a given 



The transformation |3| from emission to inertial coordinates may 
be written in a totally symmetric form without the choice of any 
emitter worldline as reference origin line. For this purpose, one 
has to consider the barycenter of the emitters as the convenient 
reference event rather than one of the emitters. This issue will 
be addressed elsewhere [26| . in connection with the symmetric 
formulation of the location problem in flat space-time. 
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configuration of the emitters. Therefore, is directly 
computable from the sole standard emission data. 

Nevertheless, if we want to obtain x from ^ we also 
need to determine the orientation e, which involves, by 
substituting ^ in the unknown x. In fact, from Eqs. 
([2]) and |7| it is clear that m a — — e a (see Fig. ^) 
and one obtains: 

mi A 777,2 A 771 3 A 777,4 = — (ei A e2 A e3 A 7774) (9) 

that taking into account ^ allows us to express Eq. ^ 
as 

e = sgn(x-m 4 ) (10) 

which by (pi) depends on x. Therefore, in order to show 
that Eq. (pF does not chase its own tail, we must be able 
to determine the orientation e at x by using a procedure 
not involving the previous knowledge of x. 




FIG. 4. For the event P, the configuration hyperplane F is 
light-like, \ 2 = 0- I n this case, the emitters remain on a 2- 
paraboloid V laying in V, and a sole emission solution P exists. 
In this 3-dimensional representation for three satellites, the 2- 
paraboloid reduces to a parabola. 



B. Analysis of the solution 

In Ref. 15,, Eq. ([3]) was obtained by separately an- 
alyzing three different cases, and gluing together their 
different solutions in a sole covariant an analytic expres- 
sion. In gluing them, the role played by the external 
element £ is essential^] 

The three cases correspond to the different causal char- 
acters of the configuration vector x- In space-time metric 
signature (— , +, +, +), one has for each case: 

(i) x time- like, % 2 < : there is a sole emission solution 
x (the other one is a reception solution) . The orientation 
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FIG. 3. For the event P, the configuration hyperplane V is 
space-like, x 2 < 0. In this case, the emitters remain on a 
2-sphere 5 laying in T, and a sole emission solution P exists. 
In this 3-dimensional representation for three satellites, the 
2-sphere reduces to a circle. 



X space-like 




FIG. 5. For the events, P and P' , the configuration hyper- 
plane r is time-like, \' 2 > 0. In this case, the emitters remain 
on a 2-hyperboloid H laying in F, and both P and P' are 
emission solutions corresponding to the same emission coor- 
dinates (t a ). In this 3-dimensional representation for three 
satellites, the 2-hyperboloid reduces to a hyperbola. 



e corresponding to the emission solution remains to be 
calculated; 

(ii) x light-like, x 2 = : there is a sole emission solu- 
tion x (the other one being degenerate). The orientation 
e corresponding to the emission solution remains to be 
calculated; 

(iii) x space-like, x 2 > : there are two emission solu- 
tions, x and x' . They only differ by their orientation e. 
The problem is how to determine the one corresponding 
to the real user. 

The above cases are illustrated in Figs. [3] [4] and [5} 
For cases (i) and (ii), the matter to determine e was 

Figure [3] shows 



For a detailed discussion about this point, see Ref. |15| . 



solved in [15] (see Sec. IV B| below), 
the emission solution for the case (i). The configuration 
hyperplane, being space-like, cuts the past light cone of 
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the solution in a 2-sphere containing the configuration of 
the emitters. Figure [4] shows case (ii), where the emitter 
configuration stays on a 2-paraboloid contained in the 
null configuration hyperplanc. 

For case (iii), e can not be determined from the sole 
emission data. Figure[5]shows a pair of emission solutions 
receiving the same emission coordinates. This 'indeter- 
mination' is known as the bifurcation problem. To solve 
it is the main subject of this paper (see Sec. |v|). 



III. THE BORDER BETWEEN THE EMISSION 
COORDINATE DOMAINS 

The emission coordinate region contains two emission 
coordinate domains (see Sec. |lV] below). The border 
between these domains is the hypersurface J, where the 
Jacobian determinant of the characteristic emission func- 
tion 9 vanishes, 



J = {x\ j e (x) = 0}. 



(11) 



We are going to obtain some related properties showing 
its interest in rclativistic positioning. 

First, let us note that, in an adequate and condensed 
form, Eq. ^ reads as 



x = 74 + V* - Ax, 



(12) 



where 
A = 



(y* ■ x) + eV A 



A = (i U ■ x? ~ vlx 2 - (13) 



As m4 = y* — ^X is a nu U vector, and vectors {y*,x} an d 
{TO4, \} generate the same 2-plane, the following relation 
holds: 



sgn (A) = sgn [(x ■ m 4 ) 2 



(14) 



and then A ^ 0, assuring consistence for the above def- 
inition of A. Consequently, one has the following result, 
made already evident by Eqs. (14 1 and (10 1. 
Proposition 1: js{ x ) — if, and only if, A = 0. 

The fact that A is non-negative says that the 2-plane 
generated by and x 1S everywhere time-like, except in 
the border J , where this plane is light-like. 

Coming back to Eq. Q , let us note that H is a simple 
bivector, that is, 



H = x A a 



for some vector a, because of i(y) * H = 0, which is 
a direct consequence of Eqs. ^ and (J6|. Therefore, 
the invariant (H, *H) vanishes, [H, *H) = 0, and the 
invariant {H,H) takes the expression (see Appendix [A|: 



On the other hand, substituting (15) into Q, is 
expressed as 

y* = ^-M-x)a-{^-a)xl (17) 



and then, Eq. (131 for A becomes 

\2 



A = (x 



2 2 

X a . 



(18) 



Consequently, A really does not depend on the choice of 
the transversal vector £ and, by comparing ( 16 ) and ( 18 1, 
the following result has been proved. 

Proposition 2: Up to sign, the quantity A defined in 



(13) is the scalar invariant (H,H) of the configuration 
bivector H: 



A = -(H,H). (19) 

Moreover, from Eq. ( |19| ), the user can determine A 
from the sole standard emission data E. Thus, taking 
into account Proposition 1, the user is able to know, from 
the sole standard set E it receives, when it is crossing the 
border J of the two emission coordinate domains. 

Furthermore, it is worth remarking that on the border 
J the location of a user may be unambiguously solved. 
There, its location is obtained from ( 12 ) by taking A = 
in Eq. Q. 

On the way, taking into account that (H, H) = 0, H 
will be a null bivector only when the invariant (H, *H) 
vanishes. Then, we have also proven the following result. 

Proposition 3: For an event x € TZ, the configura- 
tion bivector H is a null bivector if, and only if js(x) = 0. 

On the other hand, an observational method allowing 
the user to detect when it is on the border J has been 
previously studied by Coll and Pozo, who stated the fol- 
lowing result [27 1 128 ] , 

Proposition ^: The border J consists in those events 
for which any user at them can see the four emitters on 
a circle on its celestial sphere. 

This result is rather counterintuitive. When the GPS 
satellites are all near the horizon, or are all too close to- 
gether on our zenith, the error in positioning is great. It 
would seem then that the optimal conditions for a pre- 
cise location would be obtained when all the satellites are 
situated on an intermediate circle of the celestial sphere 
(say, among 30 or 60 degrees with respect to the zenith) . 
Nevertheless, Proposition 4 shows that the circle corre- 
sponds to the most degenerate distribution that a set of 
satellites may have. 

Proposition 4 also makes clear that the border J may 
be plotted from the sole observational data, a result that 
(15) was not, a priori, evident. 



IV. REGIONS AND COORDINATE DOMAINS 
IN RELATIVISTIC POSITIONING 



(H,H) 
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This section provides a geometrical background to an- 
(16) alyze the space-time regions which are relevant in rela- 
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tivistic positioning. In particular, we study the subset of 
the emission coordinate region C, where the orientation 
e is computable from the standard emission data E (the 
central region of the positioning system). 



A. Emission configuration regions C s ,Ce, and Ct 

The emission coordinate region C is constituted by 
three disjoint regions, and one can write C = C s L)CiUCt- 
They are the space-like C s , the nullCe and the time-like Ct 
emission configuration regions defined by the conditions 
X 2 < 0, x 2 — u ! an d X 2 > 0, respectively. 

This means that at every event x G C s (x £ Ci or 
x £ Ct, respectively) a user receives the signals from four 
emission events that generate a space-like (null or time- 
like, respectively) hyperplane. 

From Eqs. ([5| and ([7]), which only involve the emitter 
configuration of the standard emission data E, the user is 
able to determine the sign of x 2 - Consequently, from the 
data set E and the above definitions, the user knows in 
what configuration region, C s , Ci, or Ct, of the positioning 
system it is traveling. 



B. The central region C 



--C a UCe 



We name C c = C s U Ci the central region of the posi- 
tioning system. 

At every event x £ C c , one has u ■ x 7^ f° r anv 
future pointing time-like vector u, because % 1S n °t space- 
like in this region. Taking into account that is a 
future pointing null vector, the sign of the scalar products 
X ■ 1TI4 and u • x is the same for any future pointing time- 
like vector u, and from Eq. ( 10 ) this sign is precisely 



the orientation of the positioning system on the central 
region. More precisely, we can prove the following result: 
Proposition 5: In the central region C c , the orientation 
of a relativistic positioning system is constant, and may 
be evaluated from the sole standard emission data E: 

MxeC c , e = sgn(u-x) (20) 

where u is any future pointing time-like vector. 



Thus, from ([20]) any user in the central region is able to 
determine the orientation of the positioning system, and 
then, from Eqs. (|3|-([7|, it can obtain its own position 
x in the inertial system from the sole standard emission 
data by substituting e = sgn (u ■ x) in (J3j) . The resulting 
sign of e will be positive or negative depending on the 
time orientation of the computed vector \. 



vicinity of C c . Therefore, the border J divides the time- 
like configuration region Ct of C. In other words, the 
whole region C t cannot be recovered by a sole coordinate 
domainF] In |15j we proved the following result. 

Proposition 6: The emission coordinate region C is 
not a coordinate domain but the union of two disjoint 
coordinate domains, called the front C F , and the back 
C B emission coordinate domains, C = C F U C B . 

The front coordinate domain C F contains the central 
region C c and a proper subset Cf of the time-like con- 
figuration region, C F C Ct- This proper subset Cf is the 
part of Ct adjacent to the central region C c , so that the 
whole front domain C F ,C F = C c UCf , has, by continuity, 
constant orientation e (the same as the central region). 
However, the orientation at s 6 Cf can not be deter- 



mined from Eq. ( 20 1 , because Proposition 5 only applies 
on C c . 

The back coordinate domain C B is not a simply- 
connected domain. In fact, the region Ct is not simple- 
connected, and its leaves are constituted by pairs of 
events {x, x'} having the same emission coordinates but 
different inertial ones, defining two well differentiated re- 
gions: if x G C B , then x 1 G Cf =C t ~C B . 

To illustrate these coordinate domains, let us consider 
the simple case of a symmetric stationary positioning 
system in flat space-time. In this case, the four emit- 
ters define a regular tetrahedron, and C B is the union 
of four connected components. The common boundary 
J of the domains C F and C B is a four-leaf hypersurface 
that contains the shadows that each satellite produces on 
the signals coming from the other ones in the region C. 
The orientation e of the positioning system only changes 
across J taking different constant value on each coordi- 
nate domain. The analogous, but simpler to draw, sta- 
tionary and symmetric 3-dimensional case is illustrated 
in Fig. |6j that shows the involved configuration regions 
and coordinate domains. 



THE BIFURCATION PROBLEM. 
OBSERVATIONAL SOLUTION 



ITS 



The above results show that the standard emission 
data E are generically insufficient to locate a user of a 
positioning system in an inertial system. 

In the past, and in connection with GNSS, this problem 
was pointed out by Schmidt pQ and studied by Abel and 
Chaffee |H [S] by introducing a "bifurcation parameter" 
(equivalent to the square x 2 of the configuration vector x 
of Eq. (pi)). Afterwards, it was referred as the bifurcation 
problem [29j. 



C. Front (C F ) and back (C B ) coordinate domains 

As a consequence of Proposition 5, the Jacobian de- 
terminant does not vanish, j®(x) 7^ 0, in the immediate 



7 Remember that, for historical reasons, the coordinate domain 
of a coordinate system is an open set, but not necessarily a 
topological domain (i. e. an open and connected set). 
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FIG. 6. Quotient space S of the stationary observers of a 
stationary and symmetric relativistic positioning system in 
fiat 3-dimensional space-time. Dots 7,7' 6 S represent sta- 
tionary user worldlines, and solid dots 7 a stand for stationary 
satellites (emitters). The emission configuration regions C s , 
Ce and Ct are differently colored. The border between C B 
and Ct is J: the surface of vanishing Jacobian. Conjugate 
events, (x,x'), having the same emission coordinates but dif- 
ferent inertial ones, necessarily occur on separate parts of Ct: 
if x 6 7 C Cf , then x' £ 7' C C B . In this 3-dimensional sit- 
uation, J is just the union of the shadows that each satellite 
produces on the signals coming from the other ones in the 
region C. 



In current practical situations in present day GNSS[^] 
the bifurcation problem may be solved by hand: simply 
checking which of the two solutions satisfies an observ- 
able pertinent constraint. Thus, for example, if a user 
stays near the Earth surface the right solution is the near- 
est to the Earth radius. However, in extended GNSS or 
more general positioning systems in the Solar System, 
the bifurcation problem cannot be so easily avoided; it 
will always be present for users traveling in the time-like 
configuration region C t . 

One could think that the bifurcation problem could be 
avoided by continuity for users traveling from the central 
region C c (where they are able to calculate the position- 
ing system orientation from the standard emission data) 
to the time-like configuration region C t . But the discrete 
character of true successive location operations, and the 
fact that it suffices of a sole instant to cross the border 
3 ', also make this possibility illusory. Not only for theo- 
retical reasons, but also for future practical applications 
where the role played by Earth based coordinate systems 
could become secondary (cf. [SI HH US] ) , it is essential to 
learn to solve this important part of the location problem, 
the bifurcation problem. 



We have seen that, from the sole standard emission 
data E, the users can know the configuration region that 
they are traveling. The bifurcation problem appears 
when this configuration region is the time-like one, Ct, 
because this region is constituted by pairs of conjugate 
events, x and x' , separated by the border J ', receiving 
the same standard emission data (see Figs. [5]and|6|. 
Conjugate events belong to different (back and front) co- 
ordinate domains, of different orientation. As Eq. ^ 
shows, the knowledge of the orientation in addition to 
the data set E solves completely the bifurcation prob- 
lem. Thus, how to extend the standard emission data 
E so as to be able to determine the orientation of the 
positioning system for the user? 

We shall suppose here that, in addition to the standard 
emission data E, the users are able to observe the relative 
positions of the emitters on their celestial sphere. We 
shall denote this extended data set by E* . 

Consider an arbitrary user of unit velocity u, at the 
event x of C. With respect to this user, the null vectors 
rriA may be decomposed as 



m A = (m A ) u u + rfiA, 



(21) 



where = — u ■ m A > and Aa denote vectors of 

the proper 3-dimensional space S u of the user, rfiA £ S u 
(cf. Eq. (PHI)). 

Let us consider the unit vectors Ha = Hi a/ '(jnA)u giv- 
ing the relative directions of propagation of the signals. 
Because the vectors — Ha point to the positions of the 
emitters A, i.e., are the unit vectors along the apparent 
line of sight of the emitters A, we say that {Ha} is a set 
of observational data. It is this set of data (or any equiv- 
alent one) which, added to the standard emission data, 
is included in E* . 

By direct substitution of (21) in the expression ^ of 
e, one has 

e = sgn{*[u A (— Hi A H2 A H3 + Hi A H2 A Hi 

-Hi A n 3 A Hi + H-2 A H 3 A Hi)}} (22) 

where we have taken into account that 0^=1 771 A for 
emission vectors ttia- And because any 3- form J- in S u 
satisfies i(u) * T = *(u A J 7 ), the above expression gives 
(see Appendix [A| , 

e = sgn[(Hi,H 2 ,H 3 ) - {H 2 ,H 3 ,Hi) 

+ (r? 3 ,n 4 ,r?i) - {Hi,Hi,H 2 )\, (23) 

where the triple product is defined according to Eqs. 
dA3l and dA4b. Then, the following result holds. 



Proposition 7: The orientation e of a relativistic posi- 
tioning system is given by 



sgn[(Hi, u 2l #3) 



(24) 



with u a = H a — Hi. 



A present day GNSS allows locating only part of the interior of 
a sphere surrounding the satellite constellation. 



Thus, from the relative positions of the emitters 
on the celestial sphere of the user, we can obtain the 




FIG. 7. Observational rule to determine the orientation e at the location of a user from the relative positions of the emitters on 
its celestial sphere. In the left diagram, the visual axis is oriented towards the spherical cap that does not contain the fourth 
emitter, and then e — +1. In the middle diagram, the four emitters are on a circle of the celestial sphere of the user, and 
thus e = (Proposition 4). In the right diagram, the visual axis is oriented towards the spherical cap that contains the fourth 
emitter, and thus e = — 1. 



orientation e. For instance, if the referred emitters 1, 2, 
3 are counterclockwise aligned on a circle of the celestial 
sphere of the user and the fourth emitter is inside this 
circle, then (fii, /2 2 , fa) > 0. Then, analyzing separately 
all the possible situations we arrive to the following rule 
to obtain the orientation. 

Observational rule to determine e. For any user in 
the coordinate region C receiving the extended data set 
E* , the orientation e of the positioning system may be 
obtained as follows: 

i) consider the circle of the celestial sphere of the user 
containing the three referred emitters, a = 1, 2, 3, 

ii) turn this circle around its center in the increasing 
sense 1 —> 2 —> 3 to orient the visual axis of the user by 
the rule of the right-hand screw, 

iii) if the fourth emitter A = 4 is in the spherical cap 
pointing out by this oriented axis, then the orientation 
is e = —1, otherwise e = +1. 

By applying this observational rule, the users receiving 
the extended emission data set E* can determine the 
orientation e and, from Eq. ([3]), their position in inertial 
coordinates. 

For a better geometric comprehension of the above ob- 
servational rule, we can consider an alternative approach 
to its proof. Indeed, let us focus on the generic situation 
in which {n a } is a basis of S u . Then the solution of the 
linear system 



1,2,3) 



(25) 



is given by s = u> a L a , with the vectors L a expressed in 
terms of the know data ua as 



L a = 



e abc fib x n c 
2(ni,n 2 ,n 3 )' 



(26) 



Now, by substituting (|26|) into (|23|) we arrive to the fol- 

(27) 



lowing expression for the orientation, 

e = sgn[(l - n 4 • L)(ni,n 2 ,n 3 )] 
where L = L 1 + I? + P . 



That ua ■ L = 1 when and only when the Jacobian 
je{x) vanishes has been known since [371 HE]- From Eq. 



(27), and according to the result stated in Proposition 



4, the events of the emission coordinate region C are all 
those for which the four emitters are not aligned on a 
circle of the celestial sphere of the users at these events. 
Then, it is possible to state that the factor [n^-L— 1) in 



(271 is positive or negative if n 4 is interior or exterior, re- 



spectively, to the oriented half cone containing the three 
emitters {n a } {a = 1, 2, 3). The unit vector axis s of this 
cone is given by 



s = cos if > 0. 



(28) 



Moreover, in terms of the basis {L a } given by Eq. ( 26 1 , 
the unit axis s has the expression 



s = L cos 



(29) 



as can be directly verified. 

Therefore, a unit vector v is in the interior of the half 
cone or at its exterior if the quanti ty v ■ s is greater or 
less than cose/?, respectively, or by (29) if v ■ L > 1, or 
v ■ L < 1. Thus, by taking v = 724, from (27) one has the 
following result. 

Proposition 8: Consider the oriented half-cone con- 
taining Hi, ni and n^. If 714 is in its interior, then 
e = — sgn[(fii, r?2, ^3)]- Otherwise, e = sgn[{n\,n2,n$)]. 

From this proposition we can recover the observational 
rule by considering all the possible relative positions of 
the unit vectors ni 7 n2 7 n 3 . Fig. [7] illustrates the applica- 
tion of the rule when {fti,n2,n 3 } is a negative-oriented 
basis of S u , that is for (ni, n2,n 3 ) < 0. 

Let us remark that the relative positions of the emitters 
in the celestial sphere of a user are Lorentz invariant: by 
Lorentz transformations between users at an event, the 
diameter of the circle as well as the positions of the emit- 
ters on it may change, but their increasing sense as well 
as the interior or exterior position of the fourth emitter 
will remain unchanged. 



9 



VI. DISCUSSION AND ENDING COMMENTS 

The main result of this paper is the observational rule 
giving the orientation e of the emission coordinates for 
the user. Together with the standard emission data, it 
gives a full operational character to formula Q , allowing 
any user to obtain the coordinate transformation from 
emission coordinates to inertial ones and, in particular, 
to locate itself in inertial coordinates. 

In the central region C c , where the orientation may be 
deduced from the sole standard emission data (Proposi- 
tion 6) , both the observed and the computed orientations 
may be contrasted. 

It is worth to remark here that the sole standard emis- 
sion data allows the users to detect when they are on the 
border J separating the two coordinate domains (Propo- 
sition 1), a situation that may be also contrasted with the 
limit of the observational rule (when the four emitters are 
on a circle of the celestial sphere of each user) . In spite 
of the fact that the border does not belong to any coor- 
dinate domain, the user can also locate itself in it (taking 
A = in Eqs. p} and (pi). 



Relativistic positioning concepts have been recently 
implemented in an algorithm giving the Schwarzschild 
coordinates of the users in terms of their emission co- 
ordinates (see [23]). If the conditions of applicability of 
our rule are given (observation of the emitters), the rule 
extends the region of validity of this algorithm. 

It is important to note that, in dealing with approxi- 
mate methods, or iterative algorithms, to solve the loca- 
tion problem in weak gravitational fields, Eq. ([3]) is the 
best zero order solution to start with. 

A numerical analysis of the quantities appearing in ([3| 
has been recently implemented [30l 131] . This analysis 
provides a numerical test of the results obtained in [TS] , 
and a promising via to deal with numerical simulations 
in modeling GNSS by starting from a fully relativistic 
conception. 
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more than a matter of habit or of preference. Almost 
all our expressions have been calculated many times in 
different notations, including the index one, and the sym- 
bols in the manuscript have been chosen as the best ones 
from the above criteria. For readers for which this nota- 
tion is not usual, we indicate here the relation between 
tensor index notation and ours. 

(i) Interior or contracted "product: i(x)T denotes the 
contraction of a vector x and the first slot of a tensor 
T. For instance, if T is a covariant 2-tensor, [i(x)T] v — 

(ii) Exterior or wedge product: A. If A and B are both 
covariant (or contravariant) antisymmetric tensors, the 
wedge product A A B is the anti-symmetrized tensorial 
product A® B. For instance, for vectors x and y, one 
has 



(a: Ay) 



y»x> 



and, for a covector and a 2-form F, 

(9 A F) a p 1 = 9 a Fp 7 + OpFja + 9yF a p. 

(iii) Space-time metric: g, defined as a four dimen- 
sional Lorentzian metric, has components g^ v (det g^ v < 
0). The signature of g is taken here as (— , +, +, +). 

The scalar or inner product of two vectors x and y is 
denoted as x ■ y = g(x, y) — g iLV x ii y v — x^y^ (in particu- 
lar, x 2 = x ■ x), and it is naturally extended to bivectors, 
X and Y, according to (X,Y) = hX^Y^, Indices are 
raised or lowered by using the metric g. 

(iv) Metric volume element: n, given by 



Va/3j6 



-\J —det g e a 



where e a pjs stands for the Levi-Civita permutation sym- 
bol, e i23 = 1- 

(v) Hodge dual operator: *. Let x be a vector, H a 
2-form, and T an 3-form. Their associated Hodge du- 
als, the 3-form *x, the 2-form *H, and the 1-form *F, 
are respectively given as (*x) a/ g 7 = rja^^x' 1 , (*H) a p = 
IVafs^H^, and (*F) a = ^-q^sF 13 '' 5 . 

If x, y, z, w are space-time vectors, one has 

[*(x A y)] aP = ~T)ai3iiu(x A yY v = r] ,p^ u x fi y 1 ' ', 



Appendix A: Notation 

The notation of this paper has not been chosen for 
academic reasons but for practical ones. This notation 
allows, generically, more compact and shorter expressions 
(in occupied space and in expended time) than index no- 
tation, improving a best understanding of the formula, 
but overall, it suggests more compact and shorter calcula- 
tions. In this subject, where expressions and calculations 
are determined to become more and more complicated, 
the choice of appropriate symbols from the beginning is 



[*(x Ay A z)] a = r\ afjl& x^ y 1 z\ 

and 

* (x Ay A z Aw) = rj a p 1 sx a y^ z lr w 5 . 

(vi) Invariants associated with a 2-form H. A space- 
time 2-form H has associated two independent invariants, 
(H, H) and (H, *H), which are given as: 

(H, H) = \h^H^, (H, *H) = ^H^(*H)^. 
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(vii) Relative splitting. For an arbitrary user of a rel- 
ativistic positioning system (space-time observer of unit 
4- velocity u, u 2 = — 1), any space-time vector m may be 
written as: 

m = m u u + m (Al) 

where m u — — m ■ u and m € S u are the time-like and 
space-like components, respectively, of m relative to u, 
rh ■ u = 0. 

(viii) Induced volume on S u . The 3-dimensional Eu- 
clidean space orthogonal to u, S u , has induced volume 
element, r] u , given by rj u = -i(u)r], that is, {r] u )^s = 
—u a r] a j3 1 8- The Hodge dual operator with respect r\ u is 



denoted as * u . 

(ix) Cross and triple products in S u . Vectors in S u are 
denoted with an arrow above them. Thus, for vectors 
a, b € S u , the vector or cross product is expressed as 

a x b = *(u A a A b) = * u (a A b). (A2) 

The scalar triple product is then given by 

(a x b) ■ c = (a, b, c) = * u (a A b A c), (A3) 

or, equivalently, 

{a,b,c)u = *{aAbAc). (A4) 
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